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Abstract 


A critical issue with oxidation protected carbon/carbon composites used for spacecraft thermal 
protection is the formation of coating pinholes. In laboratory experiments, artificial pinholes were 
drilled through SiC-coatings on a carbon/carbon material and the material was oxidized at 600, 
1000, and 1400°C at reduced pressures of air. The attack of the carbon/carbon was quantified by 
both weight loss and a novel cross-sectioning technique. A two-zone .diffusion control model was 
adapted to analyze this problem. Agreement of the model with experiment was reasonable at 
1000 and 1400°C; however results at lower temperatures show clear deviations from the theory 
suggesting that surface reaction control plays a role. 


Introduction 
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Carbon/carbon composites are lightweight with desirable mechanical properties at high 
temperatures. Of course, the reactivity of carbon in oxidizing atmospheres at high temperatures 
is the drawback. To this end, several oxidation-protection schemes have been developed. These 
include silicon carbide (SiC) with an over-coating of glass (1), SiC with various boron-containing 
compounds, and SiC with an over-coating of refractory oxides (2). The SiC/glass protection 
scheme shown in Figure 1 was relevant to the current study. 

Regardless of the specific coating scheme, pinholes in the protective coating are a potential 
problem (3). These provide a path for oxygen to attack the carbon/carbon with formation of a 
cavity below the coating and eventual consumption of the underlying carbon/carbon. The purpose 
of this paper is to report observations of these cavities and present a model for their formation. 
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Carbon oxidation has been extensively studied (4,5). At high temperatures, carbon oxidizes 
primarily to CO(g) when the supply of oxygen is limited: 

C(s) + !4 0 2 (g) = CO(g) (1) 

However air contains both 0 2 (g) and C0 2 (g), which are not thermochemically compatible with 
CO(g). For this reason, there is an additional reaction step in carbon oxidation: 

CO(g) + % 0 2 (g) = C0 2 (g) (2) 

The gas-phase diffusion-controlled oxidation of a carbon substrate through pores has been 
treated by several investigators for a variety of applications. These applications include the 
burnout of carbon between two alumina plates (6), the oxidation of TaC (7,8), oxidation of carbon 
fiber coatings in composites (9-11), and the oxidation of carbon in an alumina matrix (12). In the 
case of TaC oxidation, a porous oxide forms as an outer layer and oxidation occurs through the 
pores. In all four cases, the process is described by two regions, as shown schematically for 
pinhole oxidation of SiC-coated carbon/carbon in Figure 2. Flux equations are set up for each 
region, together with boundary conditions. Solutions to the ciffusion equation lead to the carbon 
consumption rate. The major difference between the treatments is the selection of a suitable 
geometry and coordinate system. 

In this paper we will focus on SiC-protected carbon/carbon. After about twelve flights small ( < 1 
mm) pinholes have been observed on the wing-leading-edge of the US space shuttle orbiter (3). 
This study is part of a larger program to understand the origin of these pinholes and to model the 
potential damage to the underlying carbon/carbon. 


The experiments reported here were done at reduced pressures to model the atmosphere entry 
environment. These experiments were then modeled with adaptations of the well-established 
models described above. Although the experiments are with SiC-coated carbon/carbon and 
emphasize re-entry applications, these experiments and models are quite general and apply to 
other coatings of carbon/carbon in other applications. 
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Experimental 
Carbon/Carbon Material 

The material is the same as that used for the wing-leading-edge of the shuttle (Lockheed Martin 
Vought Systems, Dallas, TX). It is shown schematically in Figure 1 and described in more detail 
in reference (13). This material consists of layers of a graphitized rayon fabric, which are 
repeatedly impregnated with a liquid carbon precursor to fill the voids and then the precursor is 
converted to carbon by pyrolysis. A diffusion conversion coating of silicon carbide is then grown 
on this substrate at high temperatures by a siliconizing pack cementation process. Next a coating 
of tetraethyl orthosilicate (TEOS) is applied via a vacuum impregnation to form silicon dioxide on 
the surface to plug surface cracks and fissures. The final coating is a “Type A" sealant which 
consists of sodium silicate and small particles and fibers of SiC. The difference in thermal 
expansion of the SiC coating and the carbon/carbon substrate leads to cracks in the SiC and the 
fluid glass from the TEOS and Type A coatings is intended to fill these cracks. Samples in the 
form of 2 cm diameter disks were used for the laboratory furnace study. Larger disks, about 7 cm 
in diameter, were used for the arc-jet study. 

Artificial pinholes were made through the glass and SiC of this material with an ultrasonic drill 
(Sonic-Mill, Albuquerque, NM). Figure 3(a) shows a top view of a carbon/carbon disk with drilled 
holes and Figure 3(b) shows a cross-section of a hole close to the SiC interface with the 
carbon/carbon. The procedure to produce this cross-section is described below. For the furnace 
specimens a grid of nine holes was drilled; for the arc-jet specimen, a grid of six holes was drilled. 


Oxidation Treatments 

For the furnace exposure, disks without pinholes were first tested at each temperature to 
determine the extent of coating loss and general outgassing under the test conditions. The 
weight loss data are shown in Table I. Two to four replicate runs were done for each condition. 
The weight losses for samples with pinholes were corrected so that the weight loss per pinhole 
reflects only the oxidation of the carbon/carbon. 
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The test specimens were heat treated in a horizontal furnace with molybdenum disilicide elements 
at a reduced pressure of air, to simulate a portion of the ent-y atmosphere. A schematic diagram 
of the apparatus is shown in Figure 4. Pressure was monitored with a capacitance manometer 
and controlled with a needle valve. Temperature was monitored by a Pt/Pt-Rh thermocouple 
adjacent to the samples. Table II lists the experimental exposures. 

After exposure, the samples were analyzed both from weight change data and microstructure. A 
simple weight change indicated the amount of carbon oxidized. Only one specimen (with nine 
holes) was run per condition and the corrected weight loss per sample is shown in Table II. The 
error on these numbers reflects only the weight loss from the blanks (Table I), where multiple 
specimens were run. 

The samples were also sectioned, ground, and examined in cross-section. This was done by first 
vacuum infiltrating and mounting the specimens in epoxy. Then the specimens were sectioned 
with a slow-speed diamond saw, taking care not to cut through the pinholes. Each section was 
then re-mounted and ground until a groove appeared corresponding to a pinhole. These were 
subsequently ground to 20 pm diamond to expose the structure around the pinhole. 

In addition, one specimen was oxidized in the arc-jet facility at NASA Johnson Space Center. 
This disk with pinholes exhibited more attack than those run in the furnace. The arc-jet exposure 
was at 1 538°C and 0.05 atm for a total time of 3.5 hr. The atmosphere is air with substantial 
amounts of atoms and ions formed in the arc. Tests were interrupted for weighing at regular 
intervals. 

In most cases, the carbon/carbon below the pinhole oxidized to form an approximate semi-circular 
cross-section, indicating that the cavity resembles a hemisphere. This geometry was incorporated 
into the model. 

Results 

A summary of the experimental data is shown in Table II. Weight losses per pinhole were 
determined by taking the total weight loss and dividing by the number of pinholes through which 
oxidation occurred. Generally this was nine, but in some instances the drilled hole did not reach 
the carbon/carbon or else the hole was clogged. Representative micrographs are shown in 
Figures 5-8 from which a number of important observations can be drawn. At 600°C, there 
appears to be no attack at this scale; however there is some weight loss. The carbon/carbon 
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substrate contains two forms of carbon-the fibers and the matrix material. It is well-known that 
the matrix material oxidizes more rapidly than the fibers. The 600°C oxidative weight loss likely 
reflects the oxidation of the matrix material, which is not be evident in the low-magnification cross- 
section. Figure 9 is a micrograph of the carbon/carbon substrate (without SiC/glass coating) 
oxidized at 600°C. This illustrates the preferential attack of the matrix material. 


At higher temperatures, there is clearly more attack and the region below the pinhole exhibits a 
bowl-like appearance. At 1600°C, substantial attack of the SiC coating is observed in addition to 
oxidation of the carbon/carbon can be seen, as shown in Figure 8. No attempt was made to 
model the 1600°C exposure, since the coating system was also attacked. 

One sample was exposed in an arc-jet (sample AJ1 in Table II). Figure 10 is a plot of weight loss 
of this specimen as a function of time. This reported total weight loss reflects weight loss from 
both the coating and oxidation of the carbon/carbon. Figures 10 and 11 and Table II indicate 
substantial attack of the carbon/carbon under these conditions. 

For the furnace exposures at 1000°C and 1400°C and the arc-jet exposure, we assumed the 
attacked region was a hemisphere and calculated the weight loss due to formation of the 
hemisphere. Thus the weight loss, W, was calculated from: 

W = 4 n/6 (2A/7t) 3/2 p = 1 .45 A M (3) 

Here A is the cross-sectional oxidized area below the pinhole and p is the density of 
carbon/carbon taken to be 1.362 gm/cc. The area was measured by the standard technique of 
counting squares on graph paper. Since the oxidized area was not a perfect hemisphere, 
measurement of an area should be more accurate than measuring a single radius or average of 
several radii. The calculated weight changes are also shown in Table II. For 1000 and 1400°C 
the weight changes calculated in this way and those determined by direct weighing show 
reasonable agreement. The measured weight changes were nearly always greater; they include 
preferential oxidation of the carbon matrix material, which may extend well into the matrix but may 
not be observed in the optical micrographs. Likewise, any offset of the cross-sections from the 
centerline can only result in a lower area. 



Model 


General Approach 

As mentioned, this oxidation process can be described by a model of countercurrent diffusion in a 
pore. Following Bernstein and Koger (6), the basic assumptions in this model are: 

1 The reaction rate is controlled by diffusion of the reactants and products within the pore. 

2. All gases act ideally. 

3. The gas phase diffusion coefficients are independent of concentration. 

4. The oxidation process is described by gaseous diffusion in two regions, corresponding to the 

diagram in Figure 2(b). yj SHOULD 77 VA7~ 7& £ 
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The approach here is similar to that used by other investigators. Fluxes, J„ are defined for each 
region according to: 


J. = - D. 




. dx) 


+ c 




/-i 


(4) 


Here D, eff is the gas phase interdiffusion coefficient of the ith component, c* is its concentration, 
and x is the distance from the carbon/gas interface. The first term in equation (4) is the standard 
Fickian diffusion term, the second term is the convection term required when additional gas 
molecules are created at the boundary between Regions I and II. Gas-phase interdiffusion 
coefficients are calculated from the standard Chapman-Enskcg relationship (14): 


A- iv, =0.001853 


1 1 ] 

1/2 , v 

M 1 

l«, + mJ 

nj 


(5) 


Here M, are the masses of the subscripted molecule (in grams/mole), T is the absolute 
temperature, P is the total pressure (in atm), o A . B is the average of the two diameters, and n is the 
collision integral determined from standard tabulations (15, 16). Table III lists the calculated 
diffusivities for the conditions of interest in this study. 


Calculation of Location of Region (l)/Region (II) Boundary 
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First consider Region (I), which has only C0 2 and 0 2 counterdiffusing through nitrogen. These 
fluxes are described by: 


J'n =-■ 


j' = - 

J COj 


A, c" 




Li ~ A j 

D c 

c g;, 

L-x, 


(6) 

( 7 ) 


The parameters c, L, and x f are defined in Figure 2(b). Note that these fluxes do not include a 
convection term, since the diffusion fluxes of 0 2 and C0 2 molecules are of equal and opposite 

magnitude. We can obtain c' COi : 


D n 


c ( o, ~ 


D(o } 


( 8 ) 


Next consider Region II, which contains both C0 2 and CO. Now we have a convection term, 
since for the diffusion flux of CO is twice that for C0 2 : 


j" — _ o j" 
u vo 


( 9 ) 


The flux of C0 2 is given by equation (4): 


J" - - D 

u co 1 u co 1 


Sc ^ 

C / L C () 1 

dx 


c J" 

L co 2 d CO, 


( 10 ) 


This can be rearranged and integrated to give: 
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(11) 
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Because c COi )) c COi and we also have an expression for c c 


co 2 


J !I =- 


c r D n 




In 


V xj j 


1 + 


C T D COl 


( 12 ) 


Equation (12) is useful because it helps us derive an expression for L/x,. At the boundary 
between Regions I and II, we have the following continuity expressions: 

= J ()i J Gi = (13) 


Upon substitution into equation (12) for J " 0j , we can get the expression for L/x,: 


= 1 + 


c" D, 


O, 


c I D( () 2 


In 


1 + - 


c"D. 


o, •'o. 


C/D, 


CO, 


(14) 
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For a system with a constant — , such as air where = 0.21 , the value for L/x, is 2 123 

c r 


Calculation of Concentration Profiles and Concentrations at Boundaries 

First consider the concentration profiles for Region (I). Equation (8) indicates that the 
concentration profile for C0 2 is linear. Since the flux of 0 2 is simply the negative of the C0 2 flux, 
the concentration profile of 0 2 is also linear. These are shown in Figure 2. 

The concentration profiles of CO and C0 2 in Region II can be calculated as follows: 


J" - - D 

°(() u co 
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-CO ( J'co 
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\ ax 


(15) 
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This can be simplified and integrated to: 

( 2c r-c'co) 
(2 c r ~ c c:o)_ 


j" = - 

u (o 


( 2c T D ro '' 
x t ) 


In 


( 16 ) 


Noting that c” a »c f ’ 0 and combining with equations (13) and (14), we get a concentration 
function: 


c " = 2c 
u co 


1 - 

f c o Do 1 
1 + 1 ° 2 


v c r ^co 2 J 



(17) 


This equation gives the reduction in concentration of CO in Region (II). According to equation 
(13) the flux of C0 2 has the opposite sign and twice the magnitude at any given point. These 
concentration profiles are shown in Figure 2(b). 

Calculation of Recession Rate 

The weight loss of carbon per unit area of pinhole is given by: 


dW t 

dt 


= K 4 o 2 = 



K P_o , Dp, 
rt(l- X/ ) 


( 18 ) 


Here W c is the weight of carbon (gr), M c is the molecular weight of carbon (gr/mole), and P Q} is the 

partial pressure of oxygen (bar). The weight loss rate per pinhole is simply this quantity times the 
area of the pinhole, nr 2 . The rate of growth of a hemisphere at the bottom of the pinhole is given 
by: 
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dW K 

~dU 


= P 


dV 

dt 


= P 


f 2ns" 


dt 


(19) 


where s is the radius of the underlying hemispherical void, V is the volume, and p is the density of 
carbon. Equation (19) can be equated to equation (18) and solved for the hemispherical radius, s: 


( 3 

1/3 

‘3r 2 M c P ()i 

D() 2 t 

v2 tc p ) 


2pRT(L- 

x f ) _ 


( 20 ) 


Functional Dependence of Cavity Growth Rate 

Table IV summarizes the dependence of cavity oxidation rate on various parameters. The 
hemispherical geometry of the cavity leads to the (time)’' 3 dependence for radius increase. Only a 
slight temperature dependence is expected for a diffusion-controlled process. The independence 
on total pressure results from the fact that the pressure terms always appear as a ratio 
P(0 2 )/P(Total). This ratio remains constant in the situation discussed. If, however, the amount of 
oxygen in the atmosphere changes, then the oxidation rate would change. A change in this ratio 
would also change L/x,, as indicated by equation (14). The dependence on pinhole geometry 
follows from equation (20). Note the significant dependence on pinhole radius. 


Comparison of Theory to Experiment 

Table II lists theoretical pinhole weight changes calculated Tom the diffusion-limited models for 
each condition. A comparison of theory to experiment indicates the effectiveness of the model. 
The major source of error in these calculations is the determination of the pinhole radius. As 
shown in Figures 3b, 5-8, and 11 the difficulties in ultrasonicaHy drilling such a small hole leads to 
non-uniform pinhole radii and variations from hole to hole. This is important since the calculated 
weight loss varies as the square of the pinhole radius, as discussed above. An average pinhole 
radius of 0.02 cm was taken for these calculations. The pinhole length tended to be more uniform 
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and was taken as 0.08 cm. Uncertainties also arise from the weight losses of the specimen 
without pinholes (Table I). 

Given these uncertainties, agreement is considered reasonable for the two higher temperature 
furnace exposures, but not at the 600°C furnace exposure. Visually, there was little attack of the 
carbon (Figure 5) and the small weight gain is due to the variation in weight loss between the 
blank and the pinhole specimen. At the lower temperatures, the shortfall of the measured weight 
changes compared to the diffusion theory predictions suggests that the process is not solely 
diffusion controlled but may be influenced by the surface reaction of carbon and the oxidant. The 
selective attack illustrated in Figure 9 further supports surface reaction control of the overall 
process. 


The weight loss for the arc-jet experiment increases with time, although it is not truly linear as the 
model predicts. The amount of carbon oxidized is substantially more than that predicted from the 
model. It is likely that the diffusion equations do not adequately describe the effects on gas 
transport in the pinhole from the external high velocity stream. 


Summary 
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The oxidation through pinholes in a SiC-protected carbon/carbon material has been investigated. 
Experimental pinholes were drilled through the SiC coating and oxidation damage was quantified 
by both weight changes and the examination of cross-sections. A 'diffusion-control model for 
carbon oxidation was developed to describe this oxidation process. Agreement between the 
experimental measurements and theoretical calculations for the 1000 and 1400°C^ exposures is 
reasonable, supporting diffusion control. At 600°C the process appears to be controlled by a 
surface reaction. 
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Table II. Exposure Conditions and oxidative attack for each pinhole 

Sample Exposure I Exposure Exposure I Average, Measured cavity I Weight Change Weight change 

Time (hr) Temperature Pressure measured weight area (mm 2 ) calculated calculated from 



AJ1 3.5 1538 0.050 -210 1 17.5 +/- 11 1-106 +/- 53 
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Table III. Diffusivities of 0 2 in nitrogen (cm 2 /sec) determined from Chapman-Enskog correlation. 


Temperature (°C) 

0.0066 bar 

uiianuy WllCiailVIl. 

0.092 bar 




600 

192 

13.8 ~~ 




1000 

358 

25.7 




1400 | 

573 

41.1 




1538 


46.3 


Table IV. Dependence of cavity growth rate on various parameters. 


Parameter 


Dependence of Cavity Weight Loss 


Dependence of Cavity Radius 


Time 


cc t 


OCt 


1/3 


Temperature 


cc T 


1/2 


-p/6 


Pressure 


Independent 


Independent 


Partial Pressure 
of Oxygen 


« Pn 


* P 


1/3 


Pinhole radius 


r 273 
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